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The Haldane gap material NENP exhibits anomalies in its Knight shift, far infrared absorption and 
field-dependent gaps, which have been explained using the staggered g-tensor that occurs due to the 
low crystal symmetry. We point out that the low-temperature susceptibility is also anomalous and 
that a consistent interpretation of all data may require consideration of the Dzyaloshinskii-Moriya 
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interaction. 



I. INTRODUCTION 



Ni(C2HsN2)2N02 (C1C>4)(NENP), is one of the best-studied quasi-one-dimensional antifcrromagnets which exhibits 
a "Haldane gap" in its excitation spectrum since the atomic spins have S=l. The inter-chain coupling, J', is estimated 
^ , to be only .0004J, where J « A8K is the intra-chain coupling and the disordered phase appears to persist down to 
zero temperature. Ignoring inter-chain couplings, the standard Hamiltonian for this system consists of Heisenberg 
CZ2 | exchange plus crystal field terms: 

£ i H = ■ S j+1 + E*{S*f + E*[(S*) 2 - (Sf) 2 }}. (1.1) 

B" 

\ The crystal field interactions split the triplet magnon excitation into 3 separate modes at energies 13.6 AT, 15.7^ and 
£ ■ 29ifi 

However, various anomalies appear in the finite field behavior of NENP. The low temperature susceptibility is much 
i larger than expected from the measured gap anisotropy^ The gap does not close at the Ising transition predicted 

to occur at a finite critical fields At low T the Knight shift (local magnetic field at a nucleus) is much larger than 
expected. 4 Production of a single magnon by far infrared absorption is observed even though this is expected to 
^ ' produce only zero wave-vector excitations and a single magnon has wave-vector near tt& 

Chiba et ali^ pointed out that the Knight shift anomaly can be explained by taking into account the staggered 
part of the gyromagnetic tensor. They observed that the local crystal structure near a magnetic Ni ion has principal 
axes which are rotated from the global crystal axes and that the local principal axes take two different orientations 
for even and odd sites along a chain. The g-tensor and also the crystal field Hamiltonian are expected to align with 
' the local crystal symmetry. This implies that the g-tensor has a staggered component so that an applied uniform 
magnetic field leads to a small effective staggered field in addition to the uniform one. Because an antiferromagnet 
responds much more strongly, at low T, to a staggered field than to a uniform one, this leads to large effects at low 
T. By considering the direction of the staggered field, this theory is successful at explaining the various satellites of 
■ the proton Knight shift associated with the various inequivalent H-atoms in the unit cell. 

Mitra and Halperin&observed that this staggered field also provides a natural explanation for the field-dependence 
of the gaps. Since this staggered field is perpendicular to the uniform field it breaks the Z 2 symmetry that would 
otherwise be present and eliminates the finite field Ising transition. Using a mean field type approximation they 
Q . attempted to fit the field-dependent gaps by the estimated staggered g-tensor* 4 Furthermore, because the field is 
staggered it halves the unit cell making wave-vectors and 7r equivalent thus explaining the far infared adsorption 
. i-H , anomaly. 

So far, no explanation has been offered, as far as we know, for the anomalously large low T susceptibility. Here 
we observe that the staggered field also provides a natural explanation for this since the measured susceptibility then 
becomes a sum of uniform and staggered susceptibilities and the latter becomes quite large (but remains finite) at low 
T. However, we find that it is not possible to consistently fit the susceptibility data in terms of a staggered g-tensor 
alone. 

We also observe that another important effect has been left out of previous explanations of these anomalies. This 
is the Dzyaloshinskii-Moriya (DM) antisymmetric exchange interaction^, 

Hdm = Y^Dj-iSj xS j+1 ). (1.2) 



The low crystal symmetry of NENP permits this interaction as well as the staggered g-tensor. A convenient way of 
treating a DM interaction is to remove it by a gauge transformation. It is possible to exactly eliminate it in favor of 
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a small symmetric exchange interaction and a small perturbation to the crystal field Hamiltonian which just slightly 
change the magnon energies. However, the combination of a DM interaction and a magnetic field is less benign. 
The gauge transformation transforms the uniform field into a combination of uniform, slowly rotating uniform and 
staggered effective fields. Thus the effective staggered field has two sources which are potentially of the same order of 
magnitude. 

We calculate the susceptibility including the staggered <?-tensor and DM interaction using a Ginzburg-Landau (GL) 
mean field approach^ We also do the calculation using a type of fermionic mean field theoryiS Either approach allows 
quite good fitting of the susceptibility data at low T . A reliable determination of these parameters will probably 
await accurate numerical results on one-dimensional chains and more accurate experiments. It is possible that still 
other effects which we continue to ignore such as the staggered crystal field interaction and inter-chain couplings are 
important. Nonetheless, we expect that our basic conclusion that the DM interaction and staggered g-tensor are of 
roughly equal importance in explaining these anomalies will remain true. 

In the next section we review the crystal symmetry of NENP. Using this plus high-T susceptibility measurements 
we estimate the uniform and staggered <?-tensor. We also derive the most general form of the DM interaction allowed 
by symmetry. We then go on to discuss the low-T susceptibilities of GL and fermion models in Sec. III. In Sec. IV 
we comment on other types of experimental data and other theoretical approaches. 



II. (/-TENSOR AND DZYALOSHINSKII-MORIYA INTERACTION IN NENP 

The ethylene-diamine molecule surrounding each magnetic Ni atom in NENP has an approximate orthorhombic 
symmetry with principal axes rotated relative to those defining the crystal space group. It is convenient to describe 
this rotation in two stages. Labeling the space group axes (a, 6, c) in the conventional way we first introduce a rotation 
matrix 1Z Z which rotates by 58° about the b-axis. This defines a co-ordinate system which we label (x,y,z). (The 
chain axis, b is identified with z.) The components of the spin operators in this co-ordinate system (S a ) are related 
to those in the crystallographic system (S al ) by: 

(2.1) 

where 1Z Z is a rotation about the z-axis by 

TZ Z = I -sin0 cos0 I , (2.2) 
V l) 

A further rotation by ±6* (9 w 10°) about the y-axis, TZy, depending on sites: 

cos 9 ± sin 9 

1 
=F sin 9 cos 9 




Tl±= | 1 ) . (2.3) 



defines the local symmetry axes around a Ni site, (£,£,7/). The + or — sign occurs for even or odd sites along a Ni 
chain. (See Fig. ^) We label the corresponding spin components S", 

S" = K±S. (2.4) 
The ordinary exchange interaction in NENP is generally assumed to be of the Heisenberg form 

Hex = J Sj ■ Sj + \. (2.5) 

3 

However, there is an important symmetry breaking in the crystal field Hamiltonian. This is expected to be diagonal 
in the S" co-ordinate system, as in Eq. ^l.ip . 

H CF = ^{^(Sj") 2 + E x [(Sf f - (Sf ) 2 ]}. (2.6) 

3 

Transforming to the S coordinate system, the crystal field Hamiltonian has a diagonal uniform part and a small 
off-diagonal staggered part. We will assume that the staggered part can be ignored in what follows. The uniform 
diagonal part could then be fit to the observed magnon gaps. This implies that E z » E x > 0. 
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FIG. 1: An illustration of the three co-ordinate systems. 1 and 2 refer to the even and odd sites. 



Another important source of anisotropy, when a magnetic field is applied, is the Lande g-tensor, g. The Zeeman 
term in the Hamiltonian is written: 



H z = f i B J2h-g j S j . (2.7) 



The g-tensor is assumed to be diagonal in the (£, £, rf) basis: 



g t 

H&n,C) = I 9n I . (2. 



g c 

The uniform and staggered g-tensors in the (x, y, z) co-ordinate system are given by: 



where 



/ 9x 





s 




s 


9y 




-( 







5z / 





l(x,y,z) = g +g 

= (T^g^.oC^r 1 , (2-9) 



g% cos 9 2 + g^ sin 9 2 

ft, | (2.10) 

g$ sin 6< 2 + g c cos 9 2 



and 



(55 — 5c) sm ^ cos ' 

I . (2.11) 

(g£ — <?f ) sin 9 cos 9 

The gyromagnetic tensor in the crystallographic coordinate system (a, b, c) can be written as 

g(a,c,6) = {K z K±) m ^ {K z ll±)- 1 . (2.12) 



A. Dzyaloshinskii-Moriya interaction 

As discussed by Dzyaloshinski: 6 and Moriya 7 , an additional exchange interaction term can appear in the Hamiltonian 
which is anti-symmetric under interchanging the two sites, the DM interaction 

H DM = Y / D j .(S j xS j+1 ). (2.13) 

i 

The possible values of the DM vectors Dj can be limited by considering crystal symmetries of NENP, which at low 
temperatures are given by the space group Pn2ia 5,1Q :. First, the compound is invariant under a translation along the 
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b (or z) by two sites. This means the DM vectors are the same among the even (or odd) links. Second, the crystal 
structure is invariant under the combined operation of one site translation along the chain (b) direction and a 180° 
rotation around b. The operation acts as S°~' c — * —5*+^, 5| — > Sj +1 . This implies that D a j and D c j are staggered, 
oc (— l) 1 while Df,,j is uniform. The other symmetry operations relate sites in one chain to sites in the others, so there 
are no further restrictions on the intra-chain DM vectors. 

A nearest neighbor DM interaction in one dimension can always be eliminated by a redefinition of the spin operators 
which varies from site to site (i.e. a gauge transformation). Let us suppose that the symmetric exchange interaction 
is SO(3) invariant. Then, choosing co-ordinates so that D oc z, we may write the combined symmetric and anti- 
symmetric exchange interactions as: 

H ex = £){[(J + iDj)SfSr +1 + h.c] + JS]S? +1 } (2.14) 

3 

We may always transform this into a purely parity-symmetric exchange interaction: 



Hex = J2Wj 2 + d*(s*s* +1 + s»s» +1 ) + js;s; +1 ], (2.15) 

by a gauge transformation: 

Sj^Sje ia *. (2.16) 
When Dj = (—1) J D, the required gauge transformation simply alternates from site to site: 

On the other hand, for a uniform Dj = D, 



(-1)^(1/2) taa-\D/J). (2.17) 
aj=a-j, (2.18) 



where 



tan _1 (D/J). (2.19) 



This gauge transformation introduces a small xxz anisotropy into the symmetric exchange interaction. Its effects on 
the crystal field Hamiltonian must also be considered. If we write this, in general, as: 

H C F = S 3 E f S 3> ( 2 - 20 ) 

3,a,b 

then the effect of the gauge transformation is: 

Ej -^Kia^Ejll^iaj), (2.21) 
where H(aj) is the rotation matrix which effects the gauge transformation of Eq. (|2.16(l : 

(cos aj sin aj \ 
-sinaj cos a, . (2.22) 
l) 

Thus the principal axes of the crystal field Hamiltonian are rotated from site to site while the eigenvalues remain the 
same. For an alternating DM interaction, this is an alternating rotation which would introduce an alternating term 
in the crystal field Hamiltonian. As discussed above, such a term is expected to already be present, before the gauge 
transformation. For a uniform DM interaction, the transformed _E-tensor in the crystal field Hamiltonian rotates 
steadily along the chain. We will assume these small effects can be ignored. 

The combination of a DM interaction and an applied field leads to more important effects. Upon performing the 
gauge transformation, the g-tensor at site j is transformed as: 

g 3 ■^g j Ti-\a j ). (2.23) 



In the case of a staggered DM interaction, this leads to an alternating term in the <?-tensor even if it was not present 
before, thus adding to the effective staggered field. A uniform DM interaction leads to a rotating effective magnetic 
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field. Both staggered and uniform DM interactions can be readily treated using field theory methods. They appear 
to be approximately as important as the staggered field in explaining the various anomalies mentioned in Sec. I. We 
show that it is possible to fit the susceptibility data quite well by taking into account the staggered and uniform DM 
interactions. 

Since the DM interaction contributes the same order of magnitude to the effective staggered field as the staggered 
gyromagnetic tenso»ii*i^, we have to combine them together. For small g s and (-y), and an arbitrary direction for 

the staggered DM vector D s = (D x ,D y ,0), the staggered field can be approximated as 

h s ~g s h + (^j)D s xg u h = Ah, (2.24) 

which is just the sum of two contributions^ 1 *^. Here we have introduced another matrix A relating the total effective 
staggered field, h s to the original laboratory field h. Note that |A| << 1. 

Thus after making the gauge transformation and discarding terms which we expect to be unimportant, the Hamil- 
tonian can be written, in the (x, y, z) (S) co-ordinate system: 

H = ]T{ JS, ■ S j+1 + E Z {S*) 2 + E X [(S*) 2 - (Sp 2 ] - p B h ■ [g»K(a ■ j)}^ - (-l) j p B h s • K{a ■ j)Sj}, (2.25) 

3 

where lZ(a ■ j) is defined in Eqs. I|2.18|l and 12.2211 . The values of J ss AAK and E z w 8K have been determined from 
fitting the magnon gaps to numerical simulations. 13 - 14 E x sa OAK is extracted by the best fit of experimental data to 
the six-spin-ring model calculation^ 



III. SUSCEPTIBILITY 



A. Mean Field Results 



In the large- s ap proximation, the Heisenberg spin chain is equivalent to a field theory, the 0(3) non- linear er-model 
(see, e.g., Ref. [ija and Ref. HrI K The Hamiltonian of this model is given by 



H=(-) j dz 



where 



1 = x 1^7- 
vg at 



(0 =1) , 



The coupling constant g and the magnon velocity take the values, at s — » oo, 



(3.1) 



(3.2) 



9 = ~, v = 2Js. (3.3) 

The original spin operators are expressed in terms of the field (f> and the spin density, I as 

Sj^i-lYstfj + ll (3.4) 

If we relax the constraint of 0(3) non-linear a- model and add a repulsive <p A interaction, which is treated perturbatively, 
a much simpler theory can be obtained. Including anisotropic terms, we then phenomenologically model the low-lying 
excitations via the following bosonic quantum field theorySiii, which we refer to as the Ginzburg-Landau model (GL 
model) : 



H 



Me 



2 n< + 2 ( ^ } + 2^ 



/i B ^2 higj k Tlkl (a ■ z)ei mn (f) m Il n - fx B ^ hjlZikja • z)pk4>k + A <ft 

iklmn ik 



(3.5) 



Here eijk is the anti-symmetric tensor with €123 = 1. We assume the gaps, normalization factors and velocity are: 



A x = 15.7K, A y = 13.6K, A 2 = 29K; 
p x = p y = 1.08, p z = 1.2; u = 120iT; A = 3.7K. 



(3.6) 
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(We use units where the spacing between neighboring Ni ions along the chains is 1.) The gaps are from neutron 
scattering experiments^ and normalization factors and velocity are from numerical simulations^ 

Tsvelik proposed a fermionic field theory model 9 for NENP. We can easily include uniform DM interaction into 
this model. But there is some problem for staggered effective field (including staggered g-tensor and staggered DM 
interaction). The staggered components of the spin operators have a very complicated representation as the product 
of three Ising order (and disorder) parameter fields. Consequently, there appears to be no simple method for treating 
a staggered field in this model. If only uniform DM interaction is taken into account, the Hamiltonian can be modified 
as 



H = I dz 



k klrnnp 



where Xk is two-component Majorana fermion field 

Xk -- 



X+,k 
X-,k 



(k = 1,2,3) , 



(3.7) 



(3.6 



the sign + (-) corresponds to the right (left) movers and \ — X T 7o- 7^ (/^ = 0, 1) are chosen as 70 = <r x , 71 = ia y . 
The advantage of this model is that the field-shifted gaps agree better with neutron scattering experiments than those 
of the bosonic model of Eq. I|3.5|l . 



B. Isotropic Susceptibility: uniform and staggered 

In this subsection, we will forget about the DM interaction and crystal field terms, and discuss the uniform and 
staggered susceptibilities of isotropic Heisenberg spin-1 chain. In the isotropic case, when an external uniform magnetic 
field is applied to the system, we assume that the uniform <?-tensor is also isotropic. We set g[is = 1 in this subsection 
only. 

The Hamiltonian for this case is given by 

H = Y / [JS r S J+1 -h-S ] ]. (3.9) 

3 

Assuming the excitations are non-interacting quasi-particles and no gap anisotropy, the uniform zero-field suscep- 
tibility per spin is then 

Xu = i < (S z f >= \ < (N + - iV_) 2 >, (3.10) 

where N± are total numbers of quasi-particles with S z = ±1. 
Since JV+ and iV_ are equal in the ground state, 

Xu = ^(<N 2 + >-<N + > 2 ). (3.11) 

Write 

N + = Y,N+, k , (3.12) 



Xu = y £(< N^ k >-< N + , k > 2 ). (3.13) 



we get 

2 

fe 

So for boson and fermion distributions, we have different susceptibilities per unit length as follows: 

_ 2 t dk e ^+^/T 

Xu,B - T J 2jT ( eV /A2+t,2fc2/ T - 1)2 ' K • ) 
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2 r dk e VA^+v^kyT 
Xu ' F = TJ 2^ {e VA*+^/T^^> (3 - 15) 

where the gap is 0.4107J 18 i 19 i 20 and the velocity is 2.5Ji2a 

We plot isotropic boson, fermion, Heisenberg spin-1 chain (transfer-matrix renormalization-group method)", and 
non-linear a- model results^ in Fig. [2J We see boson model result is consistent with Heisenberg spin-1 chain result 
below around 0.2J, fermion model below roughly 0.5J, non-linear sigma model result is the best, below roughly 1.5J. 




FIG. 2: Comparison of isotropic uniform susceptibilities of different models. 



When a staggered magnetic field is applied to the system, the Hamiltonian is 

H = ^[JSj ■ S j+1 - {-iyti ■ S 3 }. (3.16) 

3 

For the free boson case, the staggered susceptibility is 

P 2 (v/J) 2 

J ■ Xs,B = (A/J) 2 = P ' ^ ' 

where p is the wave function renormalization of bosonic field. We choose p — 1.11 to fit the low temperature results 
of Heisenberg spin-1 chain (TMRG)i^i They are plotted in Fig. [3] We see they are consistent up to around 0.1J. We 
also notice this p is consistent with the average value got from the numerical simulations of equal-time correlation 
function^ (in Ref. [f| 9i= $)'• 

P= -^(Px+Py + Pz) = 1.12. (3.18) 



C. Susceptibility of GL Model: uniform and staggered 



Susceptibility data has been published for an applied field along the crystallographic a, b or c axis, which is 
anisotropic. We assume that what is measured is d 2 F/dh 2 \o etc. where F is the free energy and h a the field 
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FIG. 3: Comparison of isotropic staggered susceptibilities of different models. 



component in the a-direction. From Eqs. H2.25|) and 1|2.24J) we get in the (x, y, z) co-ordinate system: 
1 d 2 F — 1 — i- 



x tk = 



L dhidhk 



= A EfeMm(| E / dTlZ ^(a • h)U mp {a ■ n ) < SI (r)Sl (0) >) 

lm Jlj2,np 

i r p 

+A H A mk (- E ( _1 ) Jl ~ i2 y o dTH ln (a-ji)1l mp (a- j 2 ) < S£(t)S£(0) >)]. 



Defining the reduced susceptibility \ (l) a t arbitrary momentum q: 



(3.19) 



(3.20) 



then plugging the explicit form of 1Z matrix into Eq. I|3.19|l and use the translation invariance of spin correlation 
function, x %k can be written as 



1 -^X U (a) + 6 ltZ x ll (0) 



■ A K A /fe 



1 - 5 l,z -U 



X U (n + a)+S l<z x ll (n) 



(3.21) 



We expect X lfe (0) to become small at low This follows from the fact that it must vanish exponentially in the 
limit where rotational symmetry around the z-axis is exact. In this case we expect that the groundstate has Sj. = 
and that there is a finite gap, = A y to the lowest state of non-zero S^. Thus, at low T, x zz (0) cx e _Ax / T . The 
fact that in NENP, w A y suggests that this symmetry is broken only by a small amount. This small symmetry 
breaking is presumably due to the E x term in Eq. Ijl.lfl an d the DM interaction. As pointed out in Ref. |23J, the 
£? x -term leads to a splitting of the gaps of first order in E x but a T = uniform susceptibility of second order in 
E x . This suggests that x zz (0)/x xx (0) should be of order (2/15) 2 = .018. This estimate was confirmed by an explicit 
calculation using the Ginzburg-Landau field theory, reviewed above. On the other hand, the experiment obtained a 
value for this ratio of about .3. While it is possible that this just reflects errors in this rough estimate and in the 
detailed mean field calculation^ which confirmed it, it seems more likely that another explanation is required. The 
explanation could reside in impurity effects or difficulties in separating the spin susceptibility from the diamagnetic 
contribution. However, later experiments^ at lower T suggest that the impurity contribution doesn't set in until 
considerably lower T and that the data over the temperature range T > 1.7 K may be dominated by the signal from 
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the pure system. Thus we are led to consider the possibility that this discrepancy may be intrinsic. In this case, the 
obvious candidate is to include the staggered and uniform DM contributions and the staggered g-tensor at low T. But 
at high T, both x 4fe (0) an d X zfe ( 7r ) g° to 2S lk /(3T). So the staggered contribution to the susceptibility is suppressed by 
the small factors of A 2 and can be dropped. ( Since the staggered g-tensor is proportional to the difference of g-tensor 
(.9? ~ 9()i one m ay think dropping these terms will eventually affect the staggered g-tensor much. We actually did 
the calculation by keeping these terms and found the result changes very little.) We also expect the relatively small 
crystal field Hamiltonian itself to become unimportant at high T. In this limit we have: 

8 2 F 

X = ~Qfc2 = Ms 2^9ik9ikX, (3-22) 

where \ ~ * 2/(3T) at large T. Using the chain rule, the experimental measurements of susceptibility data in the 
crystallographic co-ordinate system at high T thus give us, approximately, the following results for the g-tensor in the 
(x, y, z) co-ordinate system: 

gl cos 2 + g 2 y sin 2 = (2.23) 2 
^sin^ + ^cos 2 ^ = (2.21) 2 

gl = (2.15) 2 . (3.23) 

Setting 9 = 10°, 4> = 58°, we have from Eq. (|2~TUI) 

5e = 2.20, 5^ = 2.24, ^=2.15. (3.24) 
Thus the uniform and staggered g-tensor in the [x, y, z) coordinate system are, from Eqs. (|2.1(J|) and l|2.11|l : 

2.20 

2.24 I , (3.25) 
2.15 

/ 0.008 \ 
g s = , (3.26) 
\ 0.008 / 

and matrix A can be derived from Eq. (|2.24|) . 

For the GL model, from Eq. I|3.5|l we can see there appear two extra terms when a magnetic field is applied to the 
system: 



5H = — I dzji 



X] h i^ik^u{ot ■ z)ei mn (j> rn U n + X hfTZ ik (a ■ z)p k (f> k 

iklran ik 



(3.27) 



Now we use the free field approximation (A=0) to calculate the susceptibilities of the spin chain. Expanding </>j and 
Hi (i=x, y, z) in terms of annihilation and creation operators: 



^ J „"' {exp[-i(u! ik t - kx)]a lk + h.c.}, (3.28) 



¥ 2LuJik 



Hi = ^2 J — - ■ ■ (~iu} lk ){exp[-i(uj lk t - kx)]a ik - h.c.}, (3.29) 
h 

where L is the number of spins and 



, 2v l Luj lk 

k 



uf k = A 2 +v 2 k 2 . (3.30) 
If there is no external field, the Hamiltonian becomes 

ffo = ^w,i(a^j t + i). (3.31) 
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Looking at 8H as a small term, we use first order perturbation theory for eigenstate |n>, i.e, 

|m >< m\5H\n > 



|n>-> \n>-J2' 



E„ 



E„ 



we have the finite-T formula for susceptibility 

_ 2 v e-^" 



|< n\(SH/h)\m >\ 



E„ 



(3.32) 



Define reduced uniform and staggered susceptibilities x«, \s at momentum q (q << tt) as: 

r/3 



Xu(q) = IE/ dre«W 1 -^)< 4(r)S^ 2 (0) > 



From Eq. (|3.32() we can calculate them in the GL model as follows: 



Xu(q) = TT / (tt- 



1 



27T ujiluL. 



(1 +ni + n m ) - — + (nj - n m ) — 



[i ^ 1 7^ to) 



and 



x?(«z) 



A 2 + V 2g2 ■ 



In Eq. (|3.34() . is the bosonic occupation number 



1 



[exp(uj ik /T) - 1] ' 



and 



k' = -k — q, uji=uj ik , m = n ik , (j' i = Wik<, n! i = n ik , 



(3.33) 



(3.34) 



(3.35) 



(3.36) 



(3.37) 



When the magnetic field is applied along b-axis, the uniform DM interaction will shift the staggered susceptibility by 
momentum a, but will have no effect on uniform susceptibility. So from Eq. I|3.19(l . the susceptibility along b-axis 
(h = hb) can be written as 



X b = (g^ B ) 2 X z u z (0) + \ [(A 13 ^) 2 + (A 23 /i B ) 2 ] \xT{a) + x y s y (a)] 



(3.38) 



Similarly, we can calculate the susceptibility along the a-axis (h = ha = h[(cos58°)x — (sin58°)y]). Now the only 
effect of uniform DM interaction is shifting the uniform susceptibility by momentum a, so we have 



X a = -[(ff^BC0s58 ) 2 + 5 ^Bsin58°) 2 ][xr(a) +xr(«)] +Mb(A 3 iCos58° - A 32 sin 58°)^f (0) 



(3.39) 



Up to now, we have ignored the velocity differences. The fc-integral converges at k — > oo so that it is not necessary to 
introduce an ultra-violet cut-off. Of course a physical cut-off (the lattice spacing) exists in the spin chain but, in the 
approximation A << J, including this effect makes only small corrections. On the other hand, taking into account 
the velocity differences (according to Refll3l. v x = v y = 121JT, v z — 11AK), the integrals diverge logarithmically at 
large k. This implies stronger dependence on the details of the dispersion relation at larger k and the ultra-violet cut- 
off. However, for the small velocity difference in NENP, we find the susceptibility has very weak cut-off dependence. 
Changing the cut-off from tt to 100-7T only changes \ a by about 1%. We just simply ignore this velocity difference. 
The staggered and uniform DM vectors are free parameters which can be chosen arbitrarily. We take 



2J 



-0.008, § 



-0.02, 



2J 



0.04 



(3.40) 



to fit the experimental data. We plot GL model results and experimental data in Fig. 21 The agreement is quite good 
at low T where the field theory approximations are expected to work. 



11 



43 



□ a-axis (experimental) 
O b-axis (experimental) 

- - a-axis (theoretical) 

— b-axis (theoretical) 












2 4 6 

Temperature (K) 



8 



FIG. 4: Measured susceptibility vs the GL model prediction. 



D. Susceptibility of Fermion Model: uniform and staggered 

For the fermion model, from Eq. i|3.7|) an extra term appears when a uniform effective field is applied to the system: 



5H = 



dz 



fi B ih k gl{R lm {a ■ z) 

^mnpXn'foXp 

klmnp 



(3.41) 



We can also use Eq. (|3.32|1 to compute the finite-T uniform susceptibility of this model. Let us set D z — first, i.e., 
no uniform DM interaction. Then 1Z will be identity matrix. We can easily get 



Xf = (9^b) 2 Xu,f(0) 



where in fermion model 



2-7T U>lUj' m 



,i i n Wl«m - \k\ ■ \k'\ - A;A m , .m^'-n 

(1 - n LF - n m ^ F ) ■ — ; h (ni tF - n m>F ) 



Ifcl • Ifc'l + A ; A„ 



UJl + U)' m 



LO m — Ul 



Here Ui f is the fermionic occupation number 



riikF 



and 



Similarly, 



[exp{u ik /T) + 1] ' 
k' = -k-q, u}' i = (j ik >, n i; F = n ikt F , n' itF = n ik > <F - 
X a F = (g x » B cos58°) 2 x^f(0) + (g v »Bsm58°) 2 x y u y F (0)- 



(3.42) 

(i^l^m). 
(3.43) 

(3.44) 
(3.45) 



(3.46) 



The comparison of the uniform susceptibility of bosonic GL model (D z — 0) and that of fermion model (D z = 0) is 
shown in Fig. The fermion model results are qualitatively similar to the bosonic results for the uniform susceptibility 
although the T = value for x a is smaller by about a factor of 50% in the fermionic model. This makes the agreement 
with the experimental data considerably worse before inclusion of staggered g-tensor and DM interaction. 
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FIG. 5: The comparison of the uniform susceptibility of GL model (D z — 0) and that of fermion model (D z — 0). 



Similar to the GL model, we can also include uniform DM interaction into susceptibility calculation. But as we 
discussed in Sec. IIIA, we don't know how to treat the staggered effective field. If we include uniform DM interaction 
and simply take the staggered contribution of GL model as that of fermionic model, the total susceptibilities can be 
calculated and are plotted in Fig. The DM vectors are chosen as 

D x D v D z , 

^ = -0.005, ^=-0.03, ^7= 0.07 (3.47) 

to fit the experimental data. The agreement is roughly as good as GL model. 



IV. DISCUSSION OF OTHER EXPERIMENTAL DATA AND OTHER THEORETICAL APPROACHES 



Other experimental anomalies requiring staggered (/-tensor and DM interaction for their explanation occur in the 
Knight shift and field-dependent gaps. Ignoring these small perturbations, the fermionic model is fairly successful 
at explaining the field-dependent gaps^ The GL model is accurate at low fields but captures only the qualitative 
features at higher fields where the Zeeman energy is of order the gapi^ We find that including the staggered g-tensor 
and staggered DM interaction does not significantly improve the agreement in the case of the GL model. Including 
the effect of the uniform DM interaction on the field dependent gaps is a fairly difficult problem even in the GL 
model approximation since it requires a non-linear treatment of a slowly rotating field, and we do not attempt it 
here. As remarked above there appears to be no simple way of including the staggered field in the fermion model, 
thus precluding a calculation of its effects on field-dependent gaps in that model. The field-dependent Knight shifty 1 
presents similar calculational difficulties using GL or fermion model. 

There is another low-energy effective field theory model that has been applied to NENP. This model was proposed 
by Mitra and Halperini^ The Hamiltonian is the following 



H 



dz 



E 



V 

-IT 

2 1 



A? 
2v 



ikln i 



(4.1) 



In this case we also ignore the small velocity difference which can be shown to have a negligible effect, as in the GL 

model. This model differs from the standard GL theory of Eq. (|3.5() by the factors of m the coupling to the 

magnetic field. These factors were introduced in Ref. (0), in order to obtain field dependent gaps which are the 
essentially the same as in the fermionic model and hence agree much better with experiment. 
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FIG. 6: Measured susceptibility vs the fermion model prediction. 



In mean field approximation, the uniform susceptibility is: 



2tt 



(1 + n x + n v ) 




(n x -n y )- 



(4.2) 



and similarly for \ a ■ Note that the fc-integral does not converges in the ultraviolet, even ignoring velocity differences, 

since the integrand behaves as (\J ^ — /v\k\ at large k. Thus the result will be more sensitive to the details 

of the ultra-violet cut-off than for the other models considered above. Choosing a cut-off, \k\ < ir, and ignoring the 
staggered field gives a value for x a {T — 0) which is about twice as large as that observed experimentally. Including 
the staggered field raises the theoretical result still higher making the agreement worse. Furthermore, choosing the 
arbitrary cut-off to be 107r, increases x a {T = 0) by a factor of about 2, also worsening the agreement. 

Sicling et al^, using the Lanczos algorithm and the density matrix renormalization group technique, studied the 
field-induced gaps, for a field in the z-direction, using a model containing an alternating field and alternating as well 
as uniform crystal field terms. Independent rotation matrices were assumed for these two types of alternating terms, 
rather than assuming that both g-tensor and crystal field tensors are diagonal in the (£, C, v) co-ordinate system, as 
seems likely. The DM interaction was not included. More numerical work of this type, including the DM interaction 
(and perhaps also the staggered crystal field terms) and considering the other field direction and the susceptibilities 
is needed to determine accurate values of the staggered (/-tensor and DM interactions. 

We would like to thank M. El-Batanouny for asking a question which prompted this investigation and for very 
helpful discussions on crystal symmetry. This research was supported by NSF grant No. DMR 02-03159. 
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